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Abstract. 

We discuss a progress in calculation of Feynman integrals which has been done with help of the differential equation 
method and demonstrate the results for a class of two-point two-loop diagrams. 



The idea of the differential equation method (DEM) 
(see [1,2]): to apply the integration by parts procedure [3] 
to an internal n-point subgraph of a complicated Feynman 
diagram and later to represent new complicated diagrams, 
obtained here, as derivatives in respect of corresponding 
mass of the initial diagram. 

Thus, we have got the differential equations for the 
initial diagram. The inhomogeneous terms contain only 
more simpler diagrams. These simpler diagrams have 
more trivial topological structure and/or less number of 
loops [1] and/or ends [2]. 

Applying the procedure several times, we will able to 
represent complicated Feynman integrals (FI) and their 
derivatives (in respect of internal masses) through a set 
of quite simple well-known diagrams. Then, the results 
for the complicated FI can be obtained by integration sev- 
eral times of the known results for corresponding simple 
diagrams. 

Sometimes it is useful (see [4]) to use external 
momenta (or some their functions) but not masses as 
parameters of integration. 

The recent progress in calculation of Feynman in- 
tegrals with help of the DEM. 

1. The articles [5] and [6]: 

a) The set of two-point two-loop FI with one- and two- 
mass thresholds (see Fig. 1) has been evaluated by a com- 
bination of DEM and Veretin programs for calculation of 
first terms of FI small-moment expansion. The results are 
given below (some ot them have been known before (see 
disscussions in [5])). 

b) The set of three-point two-loop FI with one- and 
two-mass thresholds has been evaluated (the results of 
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some ot them has been known before (see [6])) by a com- 
bination of DEM and Veretin programs for calculation of 
first terms of FI small-moment expansion. The finite parts 
of the integrals can be written in terms of the generalized 
Nielsen polylogarithms (see [8]). 

2. The article [7]: 

The full set ot two-loop onshell master diagrams has 
been evaluated by DEM and Kalmykov programs for cal- 
culation of first terms of FI small-moment expansion. It 
has been observed that finite parts of all such integrals 
without subdivergences can be written in terms of three 
constants, two for the real part, £3 and 7r,C7 2 (7T,/3), and 
one for the imaginary part, 7r,i^ 2 . 

3. The articles [9]: 

The set of three-point and four-point two-loop mass- 
less FI has been evaluated. 

Here we demonstrate the results of FI are displayed 
on Fig.l. 

We introduce the notation for polylogarithmic func- 
tions: 

Lia(z) = S a -i,i(z), 

S a +i,b{z) = — Til— : at. 

alb] Jo t 

We introduce also the following two variables 

m 2 ' l + v/z/(z-4) ' 

Then 

q 2 h = ~log 2 (-z)log(l -z) -21og(-z)Li 2 ( Z ) 

+3Li 3 (z) - 6S X , 2 (z) - log ( 1 - z) & + 2Li 2 (z)] , 
g 2 I 5 = 2C 2 log(l +z) +21og(-z)Li 2 (-z) 

+ log 2 ( -z) log ( 1 + z) + 4 log( 1 + z)Li 2 (z) - 2Li 3 ( -z) 



h- 




FIGURE 1. Two-loop selfenergy diagrams. Solid lines denote propagators with the mass m; dashed lines denote massless propagators. 




-2Li 3 (z) + 2S lj2 (z 2 ) - 4S L2 (z) - 45i, 2 (-z) , 
q 2 ■ 1 1 2 = Li 3 (z) - 6 ^3 - C2 log y - ^ log 3 y - 4 logy Li 2 (y ) 

+4Li 3 (y)-3Li3(-y) + iLi 3 (-y 3 ), 

9 2 -/i 3 = -6Si, 2 (z)-21og(l -z)K 2 + Li 2 (z)], 
? 2 -/i4 - log(2-z)[log 2 (l -z) -21og(-z)log(l -z) 

-2Li 2 (z)] - ^log 3 (l -z) +log(-z)log 2 (l -z) 
-2£ 2 log(l -z) -Si, 2 (l/(1 -z) 2 )+2Si, 2 (l/(l -z)) 
+25i, 2 (- 1/(1 -z)) + ilog 3 y + log 2 y [21og(l + y 2 ) 

-31og(l -y +/)] -6^3 -Li 3 (-y 2 ) + ^Li 3 (-y 3 ) 

-6Li 3 (-y) + 21ogy [Li 2 (-y 2 ) - Li 2 (-y 3 ) + 3Li 2 (-y)] , 
q 2 -I l5 = 2Li 3 (z) - log(-z)Li 2 (z) +C 2 log(l -z) 

+^log 3 y- ^log 2 y [81og(l -y) -31og(l -y+y 2 )] 

-6£ 3 - iu 3 (-y 3 ) + 3Li 3 (-y) +8Li 3 (y) 

+ logy [Li 2 (-y 3 ) - 3Li 2 (-y) - 8Li 2 (y)] , 
q 2 -hn = -C2[log(l -z) + logy] -6^3 

- 2log(l- y+y 2 ) log 2 y 

+Li 3 (-y 3 ) -9Li 3 (-y) -21ogy[Li 2 (-y 3 )-3Li 2 (-y)], 
q 2 -h 2 5 = -21og 2 ylog(l -y) -6C 3 + 6Li 3 (y) 
-61ogyLi 2 (y), 



q 2 -/i234 = -6^3 - 12Li 3 (y) -24Li 3 (-y) + 81ogy [Li 2 (y) 
+2Li 2 (-y)] +21og 2 y [log(l -y) +21og(l +y)] . 
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